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Electronic and transport properties of a short Ni nanowire suspended between two semi-infinite 
ferromagnetic Ni leads are explored in the framework of density-functional theory. The spin- 
dependent ballistic conductance of the nanowire is calculated using a scattering-based approach 
and the Landauer-Biittiker formula. The total calculated conductance in units of Go = 2e 2 /h is 
around 1.6, in fairly good agreement with the broad peak observed around ~ 1.5 for the last con- 
ductance step in break junctions. Separating contributions from different spins, we find nearly 0.5 
Go from the majority spin s-like channel, whereas the remaining minority spin conductance of 1.1 
Go contains significant contributions from several d states, but much less than 0.5 Go from s states. 
The influence of the structural relaxation on the magnetic properties and the ballistic conductance 
of the nanowire is also studied. 

PACS numbers: 73.63.Rt, 73. 23. Ad, 73.40.Cg 



INTRODUCTION 

Modern experimental techniques involving ballistic 
conductance measurements in tip-based devices or in me- 
chanically controllable break junctions allow the study 
of contacts of ultimate dimensions consisting of just a 
single atom, or of a short monatomic wire in the con- 
striction. Investigations of metal nanocontacts during 
elongation show that the conductance decreases stepwise 
following atomic rearrangements. In noble metals (Cu, 
Ag, Au) and alkali metals (Li, Na, K) the last conduc- 
tance step before breaking, most likely corresponding to 
a monatomic contact or a nanowire, takes a value close to 
Go (Go = 2e 2 /h is the conductance quantum) which 
corresponds to the free propagation of two s-like conduc- 
tion electrons, one per spin. Transition metal nanocon- 
tacts (e.g., Ni, Co, Fe) whose atoms have a partially occu- 
pied d shell, also display conductance steps but show less 
clear evidence for conductance quantization. Some con- 
ductance histograms for Ni |2j,|3j and Fe nanocontacts 
exhibit a well-defined first peak near 1 Go; some others 
show a striking fractional peak at 0.5 Go in Ni 0,0 and 
Co.Q Generally, however, low temperature break junc- 
tion measurements on Ni, Co, and Fe nanocontacts find a 
broadly peaked conductance histogram well above 1 Go. 
Untiedt et al. [|| reported a value at roughly 2 Go (in Fe) 
and 1.3 Go (in Co and Ni) while Bakker et al.f^ found a 
conductance peak at 1.6 Go in Ni. 

A variety of methods have been developed to address 
theoretically the problem of electron transport in atomic 

scale conductors.il II 113, IIH E3, EE E3, EE EE IH O 

Choi and Ihm[3 

in particular proposed a complex wave- 
vector band method particularly suitable for implemen- 
tation within standard plane- wave ab-initio density func- 
tional pseudopotential approaches. Owing to their ex- 



ponential decay, complex wave-vector Bloch states pro- 
vide the spatially localized states crucially involved in the 
lead-nanocontact-lead wave-function matching. For ap- 
plication to transition metal nanocontacts, we recently 
presented 01 an extension of Choi and Ihm's method 
to ultrasoft pseudo-potentials, making the treatment of 
very localized d states accessible with a plane-wave ba- 
sis set. The method was applied to study the effect of 
a magnetization reversal on electron transport in infinite 
monatomic Co and Ni nanowires. Elj Interestingly, these 
calculations showed that a spin reversal can completely 
block all the d channels at the Fermi energy, causing a 
large "ballistic magneto-resistance" and leaving a con- 
ductance of only ~ 1 Go of the freely propagating s-like 
states of both spins. However, the nanocontact geome- 
try could not be addressed by such a tipless calculation, 
which omits completely the finite electron reflection at 
the contact. In this paper we use the same method to 
fully include the tips, and explore the joint effects of the 
geometrical constriction and of magnetism on the ballis- 
tic transport of a magnetic metal like Ni. 

For our purpose, we consider an idealized system con- 
sisting of a three-atom Ni monatomic nanowire placed 
between two semi-infinite Ni(OOl) bulk-like "leads". Ac- 
tually, break junction experiments showed spontaneous 
formation of monatomic nanowires only for the heavy 
5d metal nanocontacts (Ir, Pt, and Au),|2(J whereas 
only transmission electron microscopy showed evidence 
of such short nanowire contacts of lighter transition met- 
als such as Co and Pd.0 Our calculation is intended as 
a case nanocontact study, rather meant to explore than 
to affirm one particular geometr y a nd its effects. We 
note that recently Bagrets et aL[21| performed similar 
calculations for the three-atom Co nanowire while Jacob 
et al.^^ studied ballistic transport in a single- atom Ni 
contacts using a cluster embedded method; [9j we shall 
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FIG. 1: Supercell used for simulating a three-atom Ni wire 
suspended between two semi-infinite Ni(OOl) leads. 



TABLE I: Structure of the three-atom Ni wire suspended 
between two Ni(OOl) leads in unrelaxed and relaxed configu- 
rations: relaxation of the first interlayer spacing in the leads 
(A12), lead-nanowire separation (D), and interatomic dis- 
tance in the nanowire (d). Some relevant theoretical data 
for the Ni(OOl) surface and for the infinite Ni nanowire are 
also provided. 





A12 (%) 


D(k) 


d(A) 


unrelaxed wire 


0.0 


1.76 


2.49 


relaxed wire 


-0.8^-2.8 


1.47 


2.19 


Ni(OOl) surface" 


-3.6 






infinite Ni wire 






2.17 



a Ref. |30| 



later compare our results with their calculations. 



CALCULATIONS: GEOMETRY, RELAXATION, 
AND MAGNETISM 

Electronic structure calculations are carried out within 
the spin density functional theory, using the standard 
plane- wave P Wscf code. [2^] The spin-polarized version of 
the generalized gradient approximation (er-GGA) in the 
form introduced by Perdew, Burke, and Ernzerhof 24] is 
used for exchange-correlation functionals. Ni atoms are 
described by ultrasoft pseudo-potentials[25j with param- 
eters given in Ref. [2^. The cutoff kinetic energies were 25 
Ry and 250 Ry for the wave functions and for the charge 
density respectively. Integration over the Brillouin Zone 
(BZ) up to the Fermi energy was performed by using spe- 
cial k pointsjl^l and a standard broadening technique [281] 
with a smearing parameter of 0.02 Ry. Periodic bound- 
ary conditions are assumed in all three directions. The 
supercell used to simulate a three-atom Ni nanowire con- 
tact is shown in Fig. A (2\/2 x 2\/2) periodicity is 
employed in the xy plane which keeps the periodically 
repeated wires sufficiently apart from one another so that 
their mutual influence can be neglected. The bulk "tips" 
are simulated by the two opposite (001) faces of a planar 
slab consisting of seven atomic (001) crystalline planes (8 
atoms per plane), quite sufficient to reproduce a bulk-like 
potential in the middle of the slab. Since the supercell is 
very large along [001] (the z axis), the BZ in this direc- 
tion is sampled only at k z — while in the (k x , k y ) plane 
perpendicular to the wire convergence demanded instead 
ten two-dimensional (2D) special k points. This sam- 
pling of the BZ was found to be sufficient for obtaining a 
converged self-consistent potential needed for subsequent 
transmission calculations. 

We started with an ideal unrelaxed atomic configu- 
ration where the end atoms of the wire are positioned 
at the ideal hollow sites of Ni(001) surfaces, their dis- 



tance D (see Fig. 1) to the surface planes being a /2 
(ao = 3.52 A is the lattice constant of the bulk Ni). 
The intra- wire spacing d was initially set equal to ao / \[2 
which is the nearest neighbor distance in the fee Ni. We 
then relaxed structurally this configuration, using a sim- 
plified (2 x 2) planar geometry (4 atoms per layer rep- 
resenting the bulk) and fifteen special k points in the 
2D BZ, which reduced the computational weight. Mini- 
mizing total energy, all atomic positions were allowed to 
relax, including the distance parameters D, d, and the 
length L of the supercell along the z axis. During re- 
laxation the wire was constrained to remain straight and 
no zig-zag-like configurations were considered. Some cal- 
culated quantities describing the geometry of the system 
are listed in Table 1. In particular, we found that the 
first interlayer distance in the leads is contracted with 
respect to the ideal unrelaxed value d\ 2 = ao/2 by value 
of 0.8 — 2.8%.|29j These values are, however, less than 
the corresponding one (3.6% d\ 2 ) calculated for the un- 
perturbed, contact free Ni(001) surface [3(1] reflecting the 
influence of the nanowire bridging two surfaces. As ex- 
pected, the most significant geometry change occurs in 
the positions of the nanowire atoms. Here, both the end 
distance D and the central spacing d contract signifi- 
cantly during the relaxation. The value of d (2.19 A) 
is in fact very close to the calculated zero-stress inter- 
atomic spacing in the infinite monatomic wire (2.17 A). 
The atomic positions obtained after relaxation were sub- 
sequently used to construct our final system with a larger 
(2^2 x 2^2) planar periodicity.^ that would have been 
much more expensive to optimize directly. 

LSDA calculations were carried out for this geometry, 
comprising a total of 59 Ni atoms. In Figs. [5^ and[2b we 
show the local magnetization after averaging in the xy 
plane as a function of z for unrelaxed and relaxed atomic 
configurations, respectively. The total magnetization in 
the lead region (marked by the dashed lines) was divided 
by 8 - the number of atoms per layer. The central layer 
magnetization is quite close to the theoretical [3(J and 
experimentally value 0.61 for magnetic moment of 
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FIG. 2: The xy integrated magnetization as a function of z 
for the tip suspended three-atom Ni wire in unrelaxed (a) and 
relaxed (b) configurations. Magnetic moments (in the Bohr 
magneton) at some atoms are also given. 

bulk Ni. The local magnetization grows when moving 
from bulk-like internal layers to the surface layer and 
further to the central atom of the nanowire. This reflects 
a common tendency in d-band metals (e.g., Ni, Co, and 
Fe) for the magnetic moment to be enhanced in condi- 
tions of reduced coordination, as found in systems with 
reduced system dimensionality where the d bands are 
narrower. The magnetic moments of the atoms at the 
tip surface and of those in interior layers (usi = 0-72 [Ib 
and nsz = 0.64 /j,b) compare well with the correspond- 
ing values calculated for the Ni(001) surface[3(3 (0.76 \j>b 
and 0.68 /is). Very reasonably, magnetization in the lead 
region was found to be quite insensitive to the structural 
relaxation. On the contrary, magnetic moments on the 
nanowire atoms reduce significantly during the relaxation 
as a consequence of a significant shortening of interatomic 
bonds. In particular, the magnetic moment on the cen- 
tral nanowire atom drops from 1.22 fis to 1.07 /ib- These 
values are very close to those (1.19 /is and 1.13 /is) for 
an infinite monatomic wire calculated at the correspond- 
ing interatomic distances. 

ELECTRON TRANSMISSION THROUGH THE 
NANOCONTACT 

Conductance is the main measurable physical quantity 
in nanocontacts. Owing to their small size, usually well 
below the electron mean free path, resistive processes of 
diffusive and of inelastic origin arc generally negligible 
in nanocontacts, where even at room temperature the 
conductance is ballistic. In the linear response regime 
(i.e., at very small applied voltages) the ballistic con- 
ductance of an open system consisting of the nanocon- 
tact connected to two semi-infinite leads is related to 



the total transmission function at the Fermi energy by 
the Landauer-Biittiker formula: G = (1/2)G T(E F ), 
where Go = 2e 2 /h is the conductance quantum. In 
our spin density functional picture electrons of different 
spin move independently in their different self-consistent 
potentials. The two spin channels are uncoupled and 
can be treated independently, so long as the spin-flip 
processes are neglected, and the total transmission is 
the sum of contributions from the two spin channels: 
T(E) — T^{E) + T±(E). In our transmission calculation, 
we consider the system shown in Fig.^ a s the nanocon- 
tact scattering region, ideally attached to semi-infinite 
leads on both sides. To obtain electron transmission at 
energy E we first construct the generalized Bloch states 
(consisting of propagating and evanescent states) which 
are solutions of the Kohn-Sham equation at energy E of 
hypothetical infinite lattice-periodic leads. These states 
are then used to construct the scattering states of the 
entire lead-nanocontact-lead system and to calculate the 
transmission function. Full details are given in Ref.|3. 

Our supercell geometry encompasses an infinite num- 
ber of nanowire replicas forming a square supercell lat- 
tice in the xy plane. The transmission (per supercell) 
of such infinite set of nanocontacts is given by an inte- 
gral over the 2D BZ of a transmission that depends on 

T(E) = j^ J T(£,k x )d 2 k ± , (1) 

where S is the area of the supercell in the xy plane. What 
we are interested in, however, is the transmission of a sin- 
gle isolated nanocontact; we should therefore try to at- 
tain the limit of an infinitely large supercell. One could 
approach this limit in many different ways, for example 
by calculating transmission at some (e.g., at the T 
point, = 0) or else by integrating the transmission 
over the 2D BZ as in Eq. 1. For a sufficiently large su- 
percell the 2D BZ will collapse to a point, so that the 
transmission function will not depend on kj_ and the two 
approaches will eventually give the same limiting result. 
We found, however, that upon increasing the supercell 
xy size the correct single nanocontact limiting value for 
transmission function is approached much faster in the 
second case. In particular, for a manageable supercell 
size such as (2y/2 x 2y/2) the kj_ point integration of the 
transmission function is crucial since it turns out to de- 
pend much on the kj_, as illustrated in Fig. [3] where the 
minority spin channel transmission of a three-atom Ni 
nanowire (in relaxed configuration) is presented for the 
different k^ values shown in the insets. The transmis- 
sion function varies significantly from one k^ to another 
reflecting a different coupling between the nanowire and 
the lead states. In particular transmission at the Fermi 
energy, which gives the conductance, changes quite no- 
ticeably with kj_. The fully kj_-integrated transmission 
(obtained by using 10 special kj^ points) is provided for 
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FIG. 3: Transmission versus energy for the minority spin 
channel of the (repeated) Ni nanowire in relaxed configura- 
tion, calculated at different kx . Selected kx values are shown 
on insets. The numerical value of transmission at the Fermi 
energy - representative of ballistic conductance at zero volt- 
age - is also given for eachi curve. The transmission function 
obtained by averaging over the full set of 10 special kx points, 
closest to the single nanowire conductance, is plotted as the 
dashed curve on the lowest panel. 

comparison in the lowest panel (dashed curve). Here kj_- 
averaging makes the transmission curve smoother, elim- 
inating several fluctuations and retaining only smoother 
features attributable to the isolated nanowire. We veri- 
fied in fact that this kx-integrated conductance depends 
much less on the supercell size than any of the selected 
kj_ points, and is thus the best approximation to that of 
the isolated nanowire. The importance of proper sam- 
pling of the 2D BZ for ballistic conductance calculation 
was also pointed out recently by Thygesen et aL[34J In 
conclusion we systematically integrated the transmission 
function over the 2D BZ of a (2^2 x 2\/2) xy supercell 
by averaging over 10 special kx points in the irreducible 
sector of the 2D BZ. [13 

CONNECTION WITH ELECTRONIC 
STRUCTURE 

Before discussing the overall ballistic conductance re- 
sults for the Ni nanowire we wish to discuss the elec- 
tronic structure of the supercell used to calculate the 
self-consistent potential. As we shall show, at each kx 
there is a clear correlation between sharp features in the 
transmission as a function of energy and the electronic 



structure of the supercell. Let us consider in particular 
the dispersion of states with k z . If the nanocontact were 
cut, all the states would be z-localized inside the lead 
block, forming a dense set of states that are flat in k z , 
totally devoid of dispersion. When the nanocontact is es- 
tablished, a majority of these states remain still localized 
inside the lead, but a few extended states appear, char- 
acterized by k z dispersion. These extended states are 
not necessarily easy to identify, because as a function of 
k z they (anti-)cross the dense lead localized states. Yet, 
those extended states that cross Ep are responsible for 
the nanocontact ballistic conductance. 

For illustration, we choose the M point of the 2D BZ. In 
Figs. 0] and [5] we present electron bands of the three-atom 
Ni nanowire (relaxed configuration) in the z direction for 
majority and minority spins, respectively. We classify the 
electron states in terms of the orbital angular momentum 
m around the nanowire axis. Though in our geometry the 
m is not strictly a good quantum number due to the pres- 
ence of the leads, all the electron states at the M point 
(as well as at the T point) can be separated into three 
groups which we call as m = 0, |77i| = 1, and |?7i| = 2. 
We do this separation as follows. At the M point (and 
at the f point), electron states are classified according to 
irreducible representations of the group C^ v . This group 
has four one-dimensional representations Ai, Ai', A2, 
and A2' and one two-dimensional representation A5. At 
nanowire atoms the {s— p 2 — d 32 2_ r 2}, {p x — d XZl p y ~ d yz }, 
{d x 2_ y 2}, and {d xy } atomic orbitals transform accord- 
ing to Ai, Ag, A2, and A2' representations, respectively. 
Therefore, we can assign m — to states of Ai, \m\ = 1 
to states of A5, and \m\ = 2 to states of A2 and A2'. 
Note that there will be in principle smaller contributions 
from higher atomic states such as /, g, but that is not 
very important. Similarly, in a transmission calculation 
each transmission eigenchannel (obtained by diagonaliz- 
ing the Hermitian matrix T + T, T being the transmission 
matrix [3^) is clearly associated with a well defined \m\ 
so that the total transmission is naturally decomposed 
onto contributions from m = 0, |m| = 1, and |m| = 2. 

Let us consider majority spin first (Fig. All the 
M point electron states as a function of k z are plotted 
by thin lines. Moreover, the states of m = symmetry 
are marked by circles when their squared wave-functions 
have a charge on the nanowire larger than 4%. The 
amount of charge for these states on both the first con- 
tact atom and the nanowire middle atom (central atom) 
is given by the height of vertical lines on the LDOS pan- 
els. In Ni all spin up d states are occupied so that the 
corresponding bands of \m\ = 1 and \m\ = 2 symmetry 
lie well below the Fermi energy. They are clearly visible 
in Fig 2] as the crowd of localized states amassed below 
energy —0.5 eV. These states should not participate to 
electron transport and we do not display in Fig. 0] the 
data for \m\ = 1 and \m\ = 2. Around the Fermi en- 
ergy the extended m = states through the wire have 
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FIG. 4: Supercell band structure for the majority spin of the 
relaxed three-atom Ni wire calculated along the z direction at 
the M point of the 2D BZ. Each electron state can be labeled 
by a certain |m|, where m is the projection of the angular 
momentum on the nanowire axis (see text). The circles mark 
the states of m = symmetry having the charge on the wire 
larger than some threshold value (see text) while the states 
of \m\ — 1 and |m| =2 symmetry are not presented since 
they are irrelevant for electron transport. For each m = 
state, the charge on the contact and central atoms of the wire 
is given by vertical line s on the LDOS panel. The electron 
transmission of the m = channel is also shown. 



mainly s — p z character. They give rise to a dispersive 
band as shown in Fig.0] Accordingly, the transmission of 
the m = channel is a smooth function of energy and is 
close to 1, because the s band is broad, and with its high 
kinetic energy it is very poorly reflected at the nanocon- 
tact. Below —0.75 eV, however, the s— p z states hybridize 
strongly with d 3z 2_ r 2 states, making their effective dis- 
persion much flatter, or the effective mass along z much 
bigger. As a result, in this energy region well below Ep 
the |m|=0 transmission drops well below 1, developing a 
strong energy dependence. 

The situation is considerably richer for minority spin 
(Fig. [SJ. Here the d states of the nanowire are par- 
tially occupied and contribute to transmission around 
the Fermi energy. Therefore, in Fig. [3] we also present 
the panels for |m| = 1 (b) and \m\ = 2 (c). We mark 
the states by circles if their charge on the nanowire is 
larger than 4% (for |m| = 1) and 16% (for |m| = 2). For 
the m = channel, the dispersive s — p z band ends well 
above Ep at an energy about 0.4 eV. The Fermi energy 
falls inside the energy region of strong s — p z — d 3z 2_ r 2 
hybridization. Here, for all \m\, the electron bands are 
quite narrow which leads to significant energy variations 
of the corresponding transmission function. As antici- 
pated there is a clear correlation between sharp features 
in the transmission as a function of energy and the pres- 
ence of extended states. For example, the maximum in 



m = 



I contact atom 

—I i-M ii HI 

central atom 




FIG. 5: The same as in Fig. [1] but for the minority spin 
channel. In addition to m = states (a), the corresponding 
data for the states of |m| = 1 (b) and of \m\ — 2 (c ) symmetry 
(which now contribute to the transmission around the Fermi 
energy) are also shown. 



the transmission for the m = channel just below the 
Fermi energy corresponds to a slightly dispersive band 
with mainly d 3z 2_ r 2 character. Similarly, the maxima 
in the |m| = 1 transmission at about —1.1, —0.5, —0.1, 
and 0.4 eV agree with the extended bands seen at those 
energies. 

It should be emphasized, however, that there are some 
aspects in the transmission function which can not be 



6 



understood solely from the band structure analysis. For 
instance, the energy band of m = symmetry has a small 
gap at an energy around —0.3 eV (majority spin) or 0.3 
eV (minority spin) where the transmission function re- 
mains continuous and does not drop to zero. This is be- 
cause the transmission function is really calculated for an 
infinite lead-nanocontact-lead open system (and should 
be in principle compared with the LDOS constructed via 
scattering states), whereas the LDOS shown refers to the 
nanocontact alone, closed onto itself via periodic bound- 
ary conditions in the z direction. The two will approach 
one another only when increasing conceptually the num- 
ber of bulk planes included in the supercell to infinity. 
The energy bands will then continuously change so that 
the band gap discussed above (which is an artifact of 
the periodic boundary condition along the z axis) will 
eventually disappear. We actually tested these conclu- 
sions using a simple tight-binding model. Note that the 
only quantity obtained from supercell calculations and 
needed for evaluation of transmission function is the self- 
consistent potential which turns out to converge much 
better (than the LDOS) with the length of the supercell. 

One feature which is important to underline again 
with regard to the m = minority spin channel is that 
its transmission at Ep is well below 1, despite a mas- 
sive presence of s states. This is clearly the result of 
a geometry-induced hybridization of s — p z states - by 
themselves highly conductive- with the poorly conduc- 
tive d 3z 2_ r 2 states. We expect this hybridization to be 
in fact more general than the assumed 3-atom nanowire 
geometry, and the s-channel transmission reduction to 
carry over to all cases where hybridization is allowed by 
symmetry. This observation will be useful later. 

BALLISTIC CONDUCTANCE 
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FIG. 6: Spin-dependent transmission function of the Ni 
nanowire in the fully relaxed configuration (upper panel) . The 
lower panels display the symmetry decomposition of the trans- 
mission function and of the LDOS at the mid-nanowire atom. 



We now present the (k^-integrated) transmission func- 
tion of the 3-atom suspended Ni nanowire in the relaxed 
configuration (Fig. [HJ . On the upper panel we show the 
total transmission for the majority (up) and minority 
(down) spin channels. To understand better the dif- 
ferent features in the transmission curves we separate 
on the lower panels the total transmission onto contri- 
butions from different \m\ and also show the LDOS at 
the middle nanowire atom of the corresponding symme- 
try. At arbitrary k^ (unlike the M or the T point), the 
transmission eigenchannels are not associated to a cer- 
tain \m\ but have, in general, contributions from all the 
|m|. Therefore we decompose the total transmission as 
follows. For each k_L and for each transmission eigen- 
channel j the charge inside a sphere of some radius tq 
(volume Vq) surrounding the middle nanowire atom is 
calculated as rij = J Vo \ipj(r)\ 2 dr. It can be further de- 
composed onto contributions from different Im (in our 
calculations I = 0, 1, and 2): rij = ^2i m n j,im, where 



Tijj m = f ° \Rj,i m {r)\ 2 r 2 dr and the projections Rj,i m (r) 
are obtained from Rj,i m (r) = J Y* m (fl)ijjj(r,tt)d£l. The 
transmission of each eigenchannel is then weighted by its 
charge of a certain symmetry. For example, the contri- 
bution of the j-th channel to the transmission of m = 
symmetry is given by Tj(l/rij)'Y2i n j.io- The contribu- 
tions from different k^ and from different eigenchannels 
are then properly added to produce the final value. 

We first discuss the transmission function for the ma- 
jority spin channel (spin up). Here, the transmission 
around the Fermi energy is almost 1 and it is due to 
nearly free propagation of the s — p z state (m = 0). Since 
the spin up d states of the wire are occupied, the \m\ = 1 
and the \m\ = 2 channels contribute to the total trans- 
mission only at energies lower than Ep . Moreover, at 
Ep the \m\ = 1 channel gives rise only to a small tun- 
neling current while the \m\ — 2 channel is completely 
closed. For the minority spin channel (spin down), the 
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transmission curves are qualitatively similar to majority 
spin channel but are shifted to higher energies due to the 
exchange gap, which in the Ni nanowire is ~ 1.25 eV. 
At the Fermi energy there are now down spin d states 
available for transport. The transmission of the m = 
channel at the Fermi energy is still close to 1 but shows 
a strong energy dependence and corresponds to a mainly 
d 32 2_ r 2 LDOS peak at the central nanowire atom. The 
transmission function for the |m| = 1 channel shows two 
broad maxima around —1.1 eV and 0.4 eV. They are as- 
sociated with two double degenerate d XZiVZ states of the 
three atom Ni wire. A third state is located at lower en- 
ergy (about —2.2 eV) and is not shown in Fig. These 
three states will develop - upon increasing the nanowire 
length - into a twofold degenerate band with bandwidth 
about 3 eV. We find that the Fermi energy is located be- 
tween two states of higher energy and the \m\ = 1 channel 
contributes here about 0.7 to the total transmission. The 
remaining \m\ = 2 channel corresponds to a very narrow 
twofold degenerate d xyx 2_ y 2 band of an infinite wire (the 
bandwidth is about 0.6 eV). This band, however, crosses 
the Fermi energy. Therefore, right below the Fermi en- 
ergy one has a noticeable contribution from the \m\ = 2 
channel to the total transmission while its transmission 
at the Ef is about 0.5. This channel is responsible for a 
transmission peak at the energy —0.1 eV. 

We can further investigate how the structural relax- 
ation affects electron transport. In Fig. we present the 
transmission function of the Ni nanowire in the unrelaxed 
initial configuration. Similarly to the relaxed configura- 
tion (Fig. , we also report the symmetry-decomposed 
transmission as well as the LDOS at the nanowire central 
atom. One can see that the main conclusions drawn for 
the relaxed configuration remain unchanged. The major- 
ity spin transmission at the Fermi energy is about 0.8 and 
is again dominated by the s-like channel. Separating the 
minority spin transmission at the Fermi energy, one ob- 
tains values of 1.0, 0.7, and 0.2 from m = 0, |m| = 1, and 
|m| = 2 channels, respectively. Due to larger wire-lead 
separation in the unrelaxed configuration (parameter D 
in Table. 1), the nanowire states are less coupled to the 
leads. As a result, the peaks in transmission curves get 
sharper (see, e.g., two peaks for the \m\ = 1 channel and 
the peak for the \m\ — 2 channel). This leads to the 
smaller transmission at the Fermi energy for the |m| = 2 
channel - 0.2 versus 0.5 in the relaxed configuration (see 
Fig. 0). For the \m\ = 1 channel, the longer interatomic 
bonds in the wire (parameter d in Table. 1) lead to the 
fact that the two maxima in the transmission are less 
separated in energy and move closer to the Fermi level. 
This compensates the previous effect so that the trans- 
mission at the Fermi energy for the \m\ = 1 channel does 
not practically change compared to the relaxed configu- 
ration. 

We can now translate all results into a calculated 
ballistic conductance. From the total transmission at 
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FIG. 7: The same as in Fig. but for the Ni nanowire in 
unrelaxed configuration. 



the Fermi energy we obtain a conductance G m (0.5 + 
1.1) G = 1.6 G and G » (0.4 + 0.95) G = 1.35 G 
for relaxed and unrelaxed nanowires, respectively. The 
conductance for the relaxed configuration turns out to be 
larger with respect to the unrelaxed case which is reason- 
able because of smaller lead-nanowire separation (better 
contact). The calculated conductances are in good agree- 
ment with break-junction measurements showing a quite 
broad first conductance peak centered at 1.6 Gq\j\ or 
at 1.3 Go.@ The breadth of the measured conductance 
histogram peak can be explained by observing that the 
detailed nanocontact conductance depends noticeably on 
the detailed geometry of the lead-nanocontact structure, 
as was seen for example by comparison of the relaxed 
and unrelaxed three-atom geometries. Recently, Bagrets 
et aL[2l| calculated the ballistic conductance of a three- 
atom Co nanowire suspended between two Co(001) leads. 
They considered an unrelaxed nanowire (which corre- 
sponds to our unrelaxed configuration) and reported val- 
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ues close to 0.5 Go for both majority and minority spin 
conductances. In comparison with our Ni results that 
finding is somewhat puzzling, as we obtain a much higher 
minority spin conductance due to a sizable contribution 
from the d electrons. In this respect, our results are closer 
to recent conductance calculations performed by Jacob 
et al.$£ They studied nanocontacts formed by two Ni 
atoms joining the two Ni electrodes and also found a 
non zero conductance from minority spin d states at low 
nanocontact stretching. 

We should also stress that conductances reported here 
are much smaller than the ideal one for an infinite tipless 
monatomic wire G = 3.5 Go, where in the absence of any 
reflection the conductance is just equal to the number of 
bands at the Fermi energy which is 7. The reduction from 
this ideal 3.5 conductance to the real 1.6 is due to the 
fact that d electrons are reflected by about 75%, and also 
to s — d hybridization. On the other hand, this limited 
transmission effect is not as drastic as the predicted ef- 
fect of a magnetization reversal (or "domain wall" ) built 
into an infinite nanowire. |l9) In the latter case all the d 
channels are completely blocked at the Fermi energy so 
that conductance of the nanowire was only G = 1 Go due 
to two (one per spin) s-like channels. 



CONCLUSIONS 

To summarize, we have presented DFT calculations of 
the structural and magnetic properties of the three-atom 
Ni nanowire suspended between two semi-infinite Ni(001) 
leads. A scattering-based approach was applied to study 
ballistic transport in both unrelaxed and relaxed config- 
urations of the nanowire. Our calculations show that 
in both atomic configurations the majority spin conduc- 
tance is attributed to one s-like channel and is, therefore, 
very close to 0.5 Go- On the contrary, the minority spin 
conductance has also noticeable contribution from the d 
states and is higher - 1.35 Go and 1.6 Go for unrelaxed 
and relaxed configurations, respectively. We found that 
overall conductances are quite close to recent experimen- 
tal break junction results. Moreover, the fact that con- 
ductance is sensitive to the details of underlying atomic 
structure could account for a quite broad first peak in 
conductance histograms observed in experiment. 
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